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Theoretical  studies  indicate  that  a  variety'  of  microscopic  and  macroscopic  plasma  instabilities  can 
significantly  affect  the  performance  of  the  magnetoplasmadynamic  (MPD)  thruster,  yet  experimental 
evidence  for  their  existence  is  extremely  limited.  The  objective  of  this  research  was  to  provide  experi¬ 
mental  evidence  for  the  existence  of  microinstabiliiies  in  a  lO-kW-class,  self-field  MPD  thruster.  The  lower- 
bound  of  the  range  of  thruster  operating  currents  investigated  in  this  study  was  set  by  the  requirement 
that  the  electromagnetic  thrust  component  be  much  greater  than  the  electrothermal  thrust  component. 
The  upper-bound  of  the  operating  current  was  set  by  the  initial  observation  of  high-frequency  thruster 
voltage  oscillations,  which  is  usually  associated  with  the  onset  of  excessive  thruster  erosion.  Experimental 
evidence  for  two  microinstabilities  was  obtained  in  the  near-field  plume  of  the  thruster:  1)  the  generalized 
lower  hybrid  drift  instability  and  2)  the  electron  cyclotron  drift  instability. 


Nomenclature 
B  =  magnetic  field 

D  =  thruster  exit  plane  diameter 

e  -  charge  of  an  electron 

/  =  frequency,  Hz 

/max  =  frequency  corresponding  to  maximum  growth  of 

the  generalized  lower  hybrid  drift  instability  in  the 
lab  frame 

/„  =  Nyquist  frequency 

/„  =  frequency  corresponding  to  maximum  growth  of 

the  electron  cyclotron  drift  instability -in  the  lab 
frame 

J  =  current  density 

/  =  thruster  current  level 

k  =  wave  vector,  2ir/A,  A  =  wavelength 
L„  =  characteristic  gradient  scale  length  of 
X  =  [V(€n  x)]-‘ 

M,  =  ion  mass 

m  =  frequency  averaging  number 

m  =  thruster  mass  flow  rate 

m,  =  electron  mass 

n  =  harmonic  number,  1,  2,  3,  ... 

ric  =  electron  number  density 

r,  6,  z  =  radial,  azimuthal,  and  axial  coordinate 

Ti  =  Larmor  radius  =.V,/f2 

rp  =  probe  radius 

5ii  =  power  spectrum 

T  =  temperature 

U  =  drift  velocity 

Vji  =.  triple  probe  voltage 

V,  =  therm^  velocity  =  i2KTIM)'° 

Ifj,  =  phase  velocity 

=  plasma  beta  = 
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K 

=  Boltzmann  constant 

=  electron  Debye  length 

V 

=  collision  frequency 

a 

=  electron  or  ion  cyclotron  frequency 

=  eB/M 

0) 

=  angular  frequency,  rad/s,  =27r/ 

=  lower  hybrid  frequency  =  (fijO,)''^ 

=  electron  or  ion  plasma  frequency  ^ 

(neVsoM) 

Subscripts 

e  =  electron 

i  =  ion 


Superscript 

-  property  at  maximum  growth  rate 
I.  Introduction 

HE  magnetoplasmadynamic  (MPD)  thruster  is  a  space 
propulsion  device  that  utilizes  the  /  X  5  force  to  accel¬ 
erate  a  plasma  to  exhaust  velocities  useful  for  missions  in 
space  (10-50  km/s).'  Currently,  the  MPD  thruster  is  not  effi¬ 
cient  enough  to  make  it  a  viable  option  for  missions  of  current 
interest.  Typical  maximum  efficiencies  are  in  the  20-40% 
range,  while  a  minimum  of  60%  is  required  for  the  MPD 
thruster  to  become  competitive  with  chemical  and  other  elec¬ 
tric  propulsion  devices.^  In  1988,  Choueiri  et  al.^  suggested 
that  current-driven  microscopic  plasma  instabilities  (i.e;,  mi¬ 
croturbulence)  may  play  an  important  role  in  the  inefficient 
operation  of  the  MPD  thruster.  Current-driven  microturbulence 
counters  the  objective  of  the  MPD  thruster,  which  is  to  push 
the  propellant,  by  taking  much  of  the  energy  available  for  ac¬ 
celeration  and  converting  it  into  various  unrecoverable  energy 
modes,  such  as  those  associated  with  thermal  motion,  ioniza¬ 
tion,  excitation,  and  radiation.  This  suggestion  provided  the 
impetus  for  this  study  of  microinstabilities  in  a  10-kW-class, 
self-field  MPD  thruster. 

An  enormous  body  of  literature  is  dedicated  to  the  study  of 
instabilities  in  plasmas  containing  a  current  perpendicular  to  a 
magnetic  field."'  There  exists  a  fair  amount  of  research  on  the 
effect  of  plasma  instabilities  on  other  electromagnetic  accel¬ 
eration  devices  utilizing  J  x  B  forces  (e.g.,  the  Hall  current 
accelerator’).  While  studies  of  plasma  instabilities  in  Hall  cur¬ 
rent  accelerators  and  other  plasiha  physics  apparatus  are  help¬ 
ful  to  MPD  thruster  research,  they  are  not  strictly  applicable, 
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because  the  plasma  properties  in  the  MPD  thruster  differ  sig¬ 
nificantly  from  those  of  other  devices.  The  plasma  internal  to 
the  MPD  thruster  is  relatively  dense  and  collisional  such  that 
the  electron  Hall  parameter  ranges  from  order  1  to  order  10 
and  the  plasma  ^  ranges  from  much  less  than  1  to  as  high  as 
order  10.  In  addition  to  the  sheath  regions,  large  gradients  in 
plasma  properties  are  typical  inside  the  MPD  thruster.  Current 
densities  and  induced  magnetic  fields  can  span  many  orders  of 
magnitude,  up  to  -10^  A/cm^  and  ~1  kG  in  MW-Ievel  de¬ 
vices.®’’  The  plasma  is  typically  highly  ionized,  with  charged 
particle  number  densities  in  the  10”-10‘®/cm^  range,®"’®  and 
is  relatively  isothermal,  with  the  electron  temperature  T,  and 
ion  temperature  T)  on  the  order  of  1  eV  (Refs.  7-13).  The 
unique  combination  of  plasma  properties  in  the  MPD  thruster 
necessitates  its  own  theoretical  and  experimental  investigation. 

Initial  work  concerning  the  effects  of  plasma  instabilities  on 
MPD  thruster  performance  focused  on  the  prediction  of  the 
onset  of  various  pathological  thruster  operation  modes.  Many 
instabilities  have  been  invoked  as  the  cause  of  the  so-called 
onset  current  limitation,  where  high-frequency  terminal  volt¬ 
age  fluctuations  and  excessive  erosion  are  observed  when  the 
thruster  current  exceeds  a  critical  level.  These  include  micro¬ 
scopic  instabilities  such  as  the  Buneman,  ion  acoustic,  and 
electron  acoustic  instabilities,’^’’®  as  well  as  macroscopic  insta¬ 
bilities  such  as  electrothermal  instability’®’”  and  ofters.’®"” 
Experimentally,  inconclusive  attempts  have  been  made  to  iden¬ 
tify  the  instability  associated  with  the  onset  current  limitation 
by  the  examination  of  power  spectra  of  both  the  terminal  volt¬ 
age  and  light  intensity  fiuctiiations.”  In  the  case  of  an  MPD 
thruster  with  an  applied  axial  magnetic  field,  plasma  instabil¬ 
ities  have  also  been  proposed  as  being  responsible  for  the  ro¬ 
tating  spoke  phenomenon.“’“ 

Only  recently  has  the  research  emphasis  shifted  towards  the 
investigation  of  the  effects  of  microinstabilities  on  "MPD 
thrusters'  operating  at  a  current  level  below  the  onset  current. 
Many  observations  support  the  view  that  microinstabilities  do 
significantly  affect  MPD  thruster  performance.  These  include 
ionization  fractions’  and  line  radiation’®  levels  much  larger 
than  what  is  expected  classically,  and  the  measurements  of 
enhanced  plasma  resistivity.”’’®  Current  research  efforts  in¬ 
clude  the  theoretical  identification  of  the  various  instabilities 
that  may  exist  in  the  MPD  thruster’’”"”’”"’®  and  the  calcu¬ 
lation  of  their  nonlinear  effects  on  MPD  thruster  perfor- 
mance.”’”’”"”  Presently,  the  focus  has  been  placed  on  the 
investigation  of  the  generalized  lower  hybrid  drift  instability 
(GLHDI),’'’’’®  because  it  is  suspected  that  this  microinstability 
most  strongly  affects  the  anomalous  transport  properties  (e.g., 
resistivity,  viscosity,  and  heat  conduction)  of  the  MPD  thruster 
plasma. 

Experimental  verification  of  the  existence  of  plasma  insta¬ 
bilities  in  the  MPD  thruster  plasma  is  extremely  limited.  Most 
notably,  by  measurement  of  the  linear  dispersion  relation,’  the 
plasma  in  the  plume  of  an  M'W-level  MPD  thruster  is  indeed 
unstable  to  small  density  disturbances.  In  addition,  recent  dis¬ 
persion  relation  measurements  have  provided  evidence  for  the 
existence  of  the  GLHDI  in  the  same  thruster.’® 

The  objective  of  this  work  was  to  provide  additional  ex¬ 
perimental  data  concerning  the  existence  of  microinstabilities 
in  the  self-field  MPD  thruster.  Specifically,  this  study  is  con¬ 
cerned  with  the  MPD  thraster  operating  at  a  sufficient  current 


level  such  that  the  electromagnetic  thrust  component  is  much 
greater  than  the  electrothermal  thrust  component, -but  at  a  level 
below  the  onset  current.  To  achieve  this  objective  three  sepa¬ 
rate  experiments  were  performed.  First,  the  plasma  properties 
in  a  steady-state  lO-kV/  MPD  thruster  plume  were  character¬ 
ized  using  a  triple  probe  and  Hall  probe.  In  the  second  exper¬ 
iment,  the  triple  probe  was  used  to  investigate  turbulent  fluc¬ 
tuations  of  plasma  properties,  as  a  means  to  identify  the 
presence  of  microinstabilities.  In  the  last  experiment,  measure¬ 
ments  of  the  plasma  dispersion  relation  were  used  to  provide 
further  evidence  for  the  identification  of  the  microinstabilities 
present  in  the  10-k'W  MPD  thruster  plume. 

This  article  is  organized  as  follows.  In  Sec.  11  a  brief  review 
of  the  microinstabilities  that  may  exist  in  the  plasma  of  the 
MPD  thruster  is  presented.  In  Sec.  HI  the  characterization  of 
the  10-kW  MPD  thruster  plume  is  discussed.  The  investigation 
of  the  turbulent  fluctuations  of  plasma  properties  and  the  mea¬ 
surement  of  the  plasma  dispersion  relation  are  discussed  in 
Secs.  IV  and  V,  respectively.  Finally,  conclusions  are  presented 
in  Sec.  'Vt. 

n.  Brief  Review  of  Microinstabilities  that  May 
Exist  in  the  MPD  Thruster  Piasma 

In  view  of  the  limited  experimental  evidence  for  the  exis¬ 
tence  of  microinstabilities  in  the  MPD  thruster,  a  brief  review 
of  those  that  may  theoretically  exist  in  the  MPD  thruster  is 
presented.  This  review  will  also  help  to  identify  regions  of 
frequency  space  most  likely  to  reveal  evidence  of  microinsta¬ 
bilities  and  to  facilitate  the  interpretation  of  experimental  re¬ 
sults  to  be  discussed  later. 

First,  it  is  important  to  specify  the  regime  of  MPD  thruster 
operation  considered  in  this  study,  which  can  be  characterized 
by  the  parameter  J^frk.  The  lower-bound  J^/m  is  established 
by  the  requirement  that  the  electromagnetic  component  of 
thrust  be  much  greater  than  the  electrothermal  component.  The 
upper-bound  J^/rh  is  limited  by  that  associated  with  the  onset 
phenomenon.  'With  these  bounds,  the  minimum  value  of 
is  estimated  to  be  -20  kA’  s/g  for  an  MPD  thmster  operating 
with  the  propellant  argon,  whereas  the  maximum  value  de¬ 
pends  on  many  parameters,  but  is  typically  (for  argon)  greater 
than  100  kA’  s/g  (Ref.  2).  MPD  thrusters  are  also  commonly 
operated  in  regimes  where  the  electrothermal  thrust  is  com¬ 
parable  with  the  electromagnetic  thrust.’®  In  such  a  regime  of 
operation,  other  macroscopic  plasma  instabilities,  not  consid¬ 
ered  in  this  study,  may  be  excited  in  the  MPD  thruster. 

In  the  previously  stated  operating  regime,  large  current  den¬ 
sities  in  the  MPD  thruster  plasma  constitute  a  large  source  of 
free  energy  to  drive  an  instability.  Experimental  measurements 
have  shown  that  the  ratio  of  the  electron  drift  velocity  to 
the  ion  thermal  velocity  Va  typically  ranges  from  1  to  100 
(Refs.  6-13).  On  this  basis,  and  for  the  sake  of  simplicit}', 
only  current-driven  microinstabilities  will  be  the  focus  of  this 
study.  The  number  of  microinstabilities  to  be  considered  is 
further  reduced  to  those  driven  by  a  cross-field  current,  be¬ 
cause  the  current  density  is  always  perpendicular  to  the  mag¬ 
netic  field  in  the  self-field  MPD  thruster. 

Listed  in  Table  1  are  five  well-known  microinstabilities 
driven  by  a  cross-field  current,  along  with  their  characteristics 


Table  1  Five  microinstabilities  known  to  exist  in  a  plasma  with  a  cross-field  current 


Instability 

Main  condition  for  existence 

Maximum  growth  characteristics  in  the  ion  rest  frame 

GLHDI 

Can  exist  for  <  Vg 

<0*  (Out 

k*  s  rZ^ 

ECDI,  CO*  =  nfl,  in  the  electron  frame 

CO*  ~  ua[(i/,/K)  -  1]"’ 

k*  -  nClAU,  -  Vg)"’ 

Buneman  instability,  co*  =  co^.  in  the 
electron  frame 

u,>v^ 

CO*  =«  {mJMif'^co^ 

k*  co^lU. 

Ion  acoustic  instability,  /crj^  »  1 

L  »  T„  U,  >  {.KTJMf^ 

CO*  ^  ftjpi 

k*  -  AJ’ 

Drift  cyclotron  instability 

<B*  nhi 

-  ri 

TILLEY  ETAL. 


383 


at  maximum  growth  (denoted  by  the  asterisk)  and  their  main 
condition  for  existence.  Further  details,  such  as  the  derivation 
of  these  characteristics,  can  be  obtained  in  the  references  cited 
later  in  this  article.  Each  instability  operates  near  a  natural 
frequency  of  the.  plasma,  and  each  has  its  own  effects  on  the 
plasma,  in  terms  of  anomalous  transport.  Microinstabilities  re¬ 
quiring  Ti  »  Tc  for  existence,  such  as  the  electron  acoustic 
instability,^’  have  been  eliminated  from  consideration,  along 
with  those  microinstabilities  that  have  a  growth  rate  much 
smaller  than  the  ion  residence  time  in  the  thruster,  e.g.,  the 
universal  instabilities."* 

The  Buneman  instability  operates  in  both  magnetized  and 
unmagnetized  plasmas  near  the  electron  plasma  frequency  (in 
the  electron  reference  frame)  and  is  excited  when  the  electron 
drift  Welocity  exceeds  the  electron  thermal  velocity.“'^*  Ex¬ 
cluding  the  sheath  regions,  experiments  in  the  MPD  thruster 
have  shown  that  the  electron  drift  velocity  is  only  a  small 

■  fraction  of  the  electron  thermal  velocity,  effectively  stabilizing 
the  Buneman  instability. 

in  a  magnetized  plasma,  for  7)  »  Tj,  the  ion  acoustic  in¬ 
stability  operates  in  two  regions  of  wave  number  space^^'^’: 
fe'Ls  «  1  (Ref.  37)  and  »  1.  In  the  region  of  kru,  »  1, 
where  significantly  higher  growth  rates  exist,  the  ion  acoustic 
instability  behaves  as  its  nonmagnetic  counterpart  with  <u* 

•  (Dpi  and  k*  ^  Aj',  where  cOpj  is  the  ion  plasma  frequency.  The 
general  condition  for  existence  of  the  ion  acoustic  instability 
is  that  TJT/ »  1,  which  is  generally  not  satisfied  in  the  MPD 
thmster.*^’'^  However,  the  ion  acoustic  instability  cannot  be 
eliminated  from'  consideration  altogether,  because  the  stability 
criterion  is  very  sensitive  to  the  actual  electron  and  ion  velocity 
distribution  functions.  Furthermore,  the  presence  of  an  electron 
temperature  gradient  has  also  been  shown  to  destabilize  ion 
acoustic  waves  in  isothermal  plasmas,"”  suggesting  that  the  ion 
acoustic  instability  may  operate  near  the  centerline  of  the  MPD 
thruster,  off  the  tip  of  the  cathode,  where  large  gradients 
.  exist.®''® 

The  drift  cyclotron  instability  (DCI)  (Ref.  41)  operates  at 
discrete  frequencies  near  the  harmonics  of  the  ion  cyclotron 
frequency  (nCl,-,  «  =  1,  2,  . . .).  The  existence  of  the  DCI  de¬ 
pends  on  the  gradient  of  various  plasma  properties  (which  cre¬ 
ate  the  cross-field  current).  The  gradient  scale  length  of  prop¬ 
erty  X,  x)],  must  be  sufficiently  small  to  excite  the 

instability,  but  not  too  small,  for  then  the  DCI  will  transform 
into  the  GLHDI  at  wave  numbers  near  Collisions  are  also 
known  to  transform  the  DCI  into  the  GLHDI  (Ref.  42).  In 
typical  MPD  thruster  plasmas,  the  gradient  scale  lengths  are 
generally  much  smaller  than  the  ion  Larmor  radius,  and  the 
ion  Hall  parameter,  which  is  the  ratio  of  the  ion  cyclotron 
frequency  to  the  ion-ion  collision  frequency,,  is  typically  much 
less  than  unity  (typically:  Clilvu  -  10“®).  Therefore,  the  DCI  is 
expected  to  be  fully  transformed  into  the  GLHDI  in  the  MPD 
thruster. 

.  The  GLHDI  (Refs.  34  and  35),  as  its  name  implies,  desig¬ 
nates  a  general  type, of  instability,  where  the  detailed  charac¬ 
teristics  (cu  and  k)  depend  on  many  plasma  parameters  and 

■  their  gradients,  but  has  the  general  tendency  to  operate  near 

the  lower  hybrid  frequency.  The  GLHDI  can  be  excited  by 
relatively  low-drift  velocities  and  has  been  studied  extensively 
in  several  more  simplified  forms,  such  as  the  modified  two- 
stream  instability,®*'®’’®®  the  lower  hybrid  drift  instability,"*®  the 
kinetic  cross-field  streaming  instability,"*"*  and  the  Farley - 
Buneman  instability."*®  The  magnitude  and  direction  of  k*  can 
vary  significantly  depending  on  the  plasma  conditions,  such  as 
the  presence  of  gradients,  collisions,  and  the  value  of  the  elec¬ 
tron  plasma  beta  The  direction  of  k*  can 

be  exactly  perpendicular  to  the  magnetic  field  or  be  on  the 
order  of  45  deg  between  the  magnetic  field  and  electron  drift 
velocity  vector;  while  the  magnitude  of  k*  can  be  as  large  as 
the  inverse  electron  gyroradius  (k*rij,  ^  I)  or  be  much  smaller 
{k*r^  «  1).  A  review  of  the  previous  special  cases  and  their 
application  to  the  MPD  thruster  plasma  is  presented  in  Ref. 


29.  The  GLHDI  apparently  has  no  fundamental  restrictions 
forbidding  its  existence  in  the  MPD  thruster.  Even  in  the 
highly  collisional  plasma  of  the  MPD  thruster,  the  GLHDI  can 
be  excited,  although  it  behaves  somewhat  differently  than  the 
collisionless  GLHDI  (Ref  29). 

In  the  electron  reference  frame,  the  electron  cyclotron  drift 
instability  (ECDI)  (Refs.  46  and  47)  operates  at  discrete  fre¬ 
quencies  near -the  harmonics  of  the  electron  cyclotron  fre¬ 
quency  (nO,  and  n  =  1,  2,  . . .).  The  typical  unstable  wave¬ 
length  is  much  shorter  than  the  electron  gyroradius  (k*ri^  » 
1),  and  is  directed  essentially  perpendicular  to  the  magnetic 
field,  in  the  direction  of  the  electron  drift  velocity.  The  ECDI, 
■like  the  GLHDI,  has  no  fundamental  restrictions  forbidding  its 
existence  in  the  MPD  thraster. 

In  summary,  excluding  unusual  circumstances  in  which  the 
ion  acoustic  instability  is  excited,  it  is  anticipated  that  two 
microinstabilities  are  most  likely  to  be  present  in  the  MPD 
thruster  plasma:  1)  the  GLHDI  and  2)  the  ECDI.  The  focus  of 
the  experiments  in  this  work  was  to  obtain  evidence  for  the 
existence  of  these  microinstabilities  in  a  10-kW  MPD  thruster. 

III.  Characterization  of  the  10-kW  MPD 
Thruster  Plume 

Knowledge  of  various  plasma  properties  is  essential  for  any 
experimental  investigation  of  plasma  instabilities.  With  the 
measurements  of  n,,  T„  the  azimuthal  magnetic  field  Bg,  and 
the  current  density  J,  it  is  possible  to  calculate  most  of  the 
relevant  plasma  properties  required  to  investigate  the  stability 
of  the  MPD  thruster  plasma.  In  this  section,  the  results  of  the 
plasma  property  characterization  in  the  plume  of  the  10-kW 
MPD  thruster  are  presented. 

All  experiments  were  conducted  in  the  near-field  plume  of 
a  steady-state  10-kW-class  MPD  thruster;  i.e.,  with  z/D  <  1, 

■  where  z  is  the  axial  distance  downstream  from  the  exit  plane. 
The  fact  that  these  experiments  were  not  conducted  inside  of 
the  thruster  did  not  'severely  limit  the  scope. of  this  investiga¬ 
tion;  many  studies  have  shown  that  significant  electromagnetic 
interactions,  characteristic  of  inside  of  the  MPD  thmster,  can 
occur  in  the  plume  (e.g..  Ref.  48).  Moreover,  concerning  the 
question  of  the  existence  of  microinstabilities,  it  has  been  sug¬ 
gested  that  the  plasma  in  the  plume  is  quite  similar  to  that 
inside  the  MPD  thraster. ‘'® 

The  experiments  were  conducted  in  the  Princeton  30-kW 
steady-state  thruster  facility.®'®®  The  vacuum  tank  is  1.5  m  in 
diameter,  6.4  m  long,  evacuated  by  one  1.2-m  diffusion  pump, 
and  backed  by  a  Roots  blower  and  a  mechanical  pump.  With 
this  pumping  system  the  tank  pressure  during  thruster  opera¬ 
tion  was  less  than  5  X  10“^  torr.  Probe  positioning  was 
achieved  by  a  pendulum  arm  mechanism,  which  swept  the 
probes  through  the  plume,®®  allowing  for  radial  profiles  of  the 
plasma  properties  to  be  obtained.  The  steady-state  MPD 
thruster  used  in  this  research  is  shown  schematically  in  Fig.  1. 

Graphite 


Fig.  1  Scale  schematic  of  the  steady-state  10-kW  self-field  MPD 
ffaruster. 


The  thruster  has  coaxial  geometry,  a  1%  thoriated  tungsten 
0.64-cm-diam  cathode,  a  graphite  anode  insert  with  an  exit 
diameter  of  3.5  cm  enclosed  in  a  graphite  cylinder  10  cm  in 
diameter.  The  plumbing  and  electrical  connections,  located  be¬ 
hind  the  thruster,  are  described  in  Ref.  50.  The  thruster  was 
powered  by  a  modified  50-kW  welding  supply  with  a  maxi¬ 
mum  ripple  factor  of  0.5%.  Experiments  were  performed  ap¬ 
proximately  30-60  s  after  ignition  to  allow  the  thruster  to 
equilibrate  thermally.  Typical  operating  conditions  with  the 
propellant  argon  are  currents  from  500-900  A,  mass  flow  rates 
in  the  8-13  mg/s  range,  and  terminal  voltages  of  20-30  V. 

A  triple  probe^‘  was  employed  to  measure  the  electron  tem¬ 
perature  and  electron  number  density  of  the  plasma.  The  triple 
probe  is  capable  of  making  essentially  instantaneous  measure¬ 
ments  of-r,  and  ji„  allowing  for  the  attainment  of  an  entire 
radial  profile  in  a  single  swing  of  the  pendulum  arm  through 
the  plume.  The  triple  probe  consisted  of  three  cylindrical  elec¬ 
trodes  configured  in  a  three-in-line  fashion  (the  electrode  spac¬ 
ing  was  1  mm,  the  was  0.121  mm,  and  the  electrode  length 
was  5  mm).^’  The  thin  sheath  limit  criterion,  which  applies  for 
TpIXd  >  100,  was  not  satisfied  in  the  plume  of  the  10-kW 
thruster  (typically,  TpIXn  -  10-30).  Therefore,  a  modification 
of  the  triple-probe  theory,'",  based  on  the  calculations  of  Laf- 
ramboise,“  was  utilized.  Using  this  theory,  and  were 
calculated  from  the  two  triple-probe  outputs:  1)  the  voltage 
and  2)  the  current.  In  these  calculations  it  was  assumed  that 
Ti  -  T,  and  that  the  ion  species  were  singly  ionized.  The  effect 
of  ions  drifting  perpendicular  to  the  triple  probe,  known  to 
significantly  affect  triple^probe  performance, was  accounted 
for  by  aligning  the  triple  probe  with  the  ion  flow  vector  as 
discussed  in  Refs.  10  and  49.  The  uncertainty  associated  with 


Fig.  2  Typical  radial  profiles  of  electron  temperature  in  the  near- 
field  plume  of  the  10-kW  MPD  thruster  (uncertainty  in  the  mea¬ 
surements  was  >=“10%). 


Fig.  3  Typical  radie!  profiles  of  electron  number  density  in  the 
near-field  plume  of  Che  10-kW  MPD  thruster  (uncertainty  of  the 
measurements  was  ’«60%). 


the  triple-probe  measurements  was  estimated  to  be  =’10%  for 
7)  and  =’60%  for  due  to  untoown  ion  temperature,  un¬ 
known  species  composition,  and  because  of  the  existence  of  a 
nonzero  electron  drift  velocity  perpendicular  to  the  probe.'" 

In  Figs,  2  and  3,  three  radial  profiles  of  electron  temperature 
and  electron  number  density  in  the  near- field  plume  of  the  10- 
kW  thruster  are  shown.  The  T,  profiles,  -shown  in  Fig.  2,  are 
very  similar  in  shape  and  magnitude  to  spectroscopic  mea¬ 
surements  obtained  by  Myers"  at  the  exit  plane  of  the  same 
thruster.  Both  indicate  a  relatively  flat  profile  of  except  near 
the  centerline.  The  magnitudes  of  the  number  density  mea¬ 
surements,  shown  in  Fig.  3',  were  also  consistent  with  spectro¬ 
scopic  measurements  obtained  at  the  exit  plane  by  Myers."  The 
effect  of  mass  flow  rate  on  T,  is  small,  while  n,  is  observed  to 
increase  slightly  with  increasing  mass  flow  rate,  as  observed 
in  other  MPD  thrusters.”  Note  that  the  triple-probe  output  was 
not  filtered,  indicating  that  the  plasma  was  not  grossly  turbu¬ 
lent  at  these  operating  conditions. 

Measurements  of  file  azimuthal  magnetic  field  and  current 
density  in  the  plume  of  the  10-kW  thruster  were  obtained  by 
the  use  of  a  Hall-effect  probe.*^  It  consisted  of  two  Hall  gen¬ 
erators,  displaced  axially  such  that  a  radial  profile  of  Be  and 
dBe/dz  could  be  obtained."*"  The  axial  and  radial  current  den¬ 
sities,  used  to  derive  the  electron  drift  velocity  with  respect  to 
the  ions  =  Jln^e),  were  then  calculated  via  Maxwell’s  equa¬ 
tions.  Figure  4  shows  typical  measurements  of  the  azimuthal 
magnetic  field  at  an  axial  location  of  1.7  cm  downstream  of 
the  exit  plane,  at  a  mass  flow  rate  of  10  mg/s,  as  the  thruster 
current  is  varied.  It  was  observed  that  the'magnetic  fields  were 
strongly  affected  by  the  thruster  current  level,  while  being  rel¬ 
atively  insensitive  to  inass  flow  rate  variations.  The  small 
hump  occurring  near  r  =  -t-5  cm,  caused  by  mechanical  factors, 
was  eliminated  in  the  experiments  to  be  discussed  in  Secs,  IV 
and  V.  . 

The  largest  source  of  error  in  the  magnetic  field  measure¬ 
ments  is  due  to  the  uncertainty  of  the  location  of  the  probe 
(Ar  =  ±0.5  steps  =  1.3  mm).  Depending  on  the  region  of  the 
plume,  the  uncertainty  in  the  magnetic  field  measurement  was 
typically  10-25%  and  no  less  than  3%,  which  was  the  cali¬ 
bration  uncertainty.  For  radii  less  than  0.75  cm,  it  is  also  ex¬ 
pected  that  the  probe  will  significantly  perturb  the  magnetic 
field  by  inhibiting  the  current  inducing  Bg  (Ref.  49).  Therefore, 
in  the  experiments  to  be  discussed  later,  only  measurements 
from  the  r  >  0.75-cm  region  of  the  plume  and  from  reasonably 
symmetric  profiles  were  considered  reliable.  Typically,  in  the 
region  r  >  0.75,  axial  current  densities  ranged  from  1  to  10  A/ 
cm^  while  maximum  radial  current  densities  were  in  the  1  -4 
A/cm"  range.  Typical  uncertainties  in  these  calculations,  be¬ 
cause  of  the  uncertainty  in  the  position  of  the  probe,  were  25- 
50%. 


Radius  (cm) 

Fig.  4  Typical  radial  profiles  of  the  asmuithal  magnetic  field  in 
the  aear-fleld  plume  of  the  lO-kW  MPD  thruster. 
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IV.  Investigation  of  Ibrbulent  Power  Spectra 

The  rationale  behind  this  experiment  stems  from  the  as¬ 
sumption  that  each  microinstability  listed  in  Table  1  operates 
near  a  different  natural  frequency  of  the  plasma.  By  measuring 
the  power  spectrum  of  naturally  occurring  plasma  fluctuations 
and  comparing  it  with  measured  natural  frequencies  (j^,  /ui, 
andj^),  it  is  possible  to  obtain  evidence  indicating  that 
microinstabilities  are  present  in  the  MPD  thruster  plasma.  Pre¬ 
vious  researchers,  in  an  attempt  to  identify  plasma  instabilities, 
have  investigated  the  power  spectra  of  both  the  MPD  thmster’s 
terminal  voltage  and  light  intensity  measurements  from  the 
plume,  without  the  concurrent  measurement  of  the  natural 
plasma  frequencies.^  Moreover,  the  terminal  voltage  and  light 
intensity  measurements  in  that  work  were  integrated  over 
space,  which  make  them  exfxemely  difficult  to  use  for  the  iden¬ 
tification  of  plasma  instabilities.  This  experiment  differs  from 
past  research  in  that  the  frequency  content  of  the  local  plasma 
property  fluctuations  are  measured  and  compared  with  mea¬ 
surements  of  local  natural  frequencies,  allowing  for  a  much 
clearer  interpretation  of  the  spectra. 

Three  main  assumptions  form  the  basis  for  this  approach. 
The  first  assumption,  consistent  with  the  reasoning  in  Sec.  E, 
is  that  Table  1  represents  an  extensive  list  of  all  microinsta¬ 
bilities  that  may  be  excited  in  the  MPD  thruster  (operating  at 
sufficiently  high  JfM).  Second,  the  characteristics  of  each  mi¬ 
croinstability  listed  in  Table  1  are  approximately  correct  for 
the  plasma  of  the  10-kW  thruster  plume,  where  jS,  is  generally 
much  larger  than  one.  This  assumption  is  reasonable,  consid¬ 
ering  that  theoretical  studies  have  shown  that  the  characteris¬ 
tics  of  the  GLHDI  (Refs.  29  and  30)  and  the  ECDI  (Ref.  55) 
do  not  change  much  as  /3,  approaches  unity.  While  it  is  true 
that  01  does  affect  /c*,  co*,  and  the  stability  of  the  microinsta- 
bilitieS  listed  in  Table  1,  only  order-of-magnitude  estimates  of 
and  a>*  are  required  to  interpret  the  data.  This  is  the  case 
because  the  following  condition  applies  to  the  plasma  in  the 
plume  of  the  10-kW  MPD  thruster: «/lh  «/«  «/pi  « 
The  third  assumption  is  that  the  location  of  the  peak  in  the 
power  spectra  is  approximately  equal  to  the  frequency  corre¬ 
sponding  to  the  maximum  linear  growth  rate.  Such  an  as¬ 
sumption  is  reasonable  because  it  is  commonly  observed  that 
saturated-state  plasma  fluctuations  still  retain  much  of  their 
linear  characteristics.^^ 

The  triple  probe  was  the  primary  means  by  which  plasma 
property  fluctuations  were  measured.  The  fluctuations  of  'Kn 
reflect  essentially  fluctuations  of  Tj  and  the  electric  field  com¬ 
ponent  along  the  electrodes,  whereas  fluctuations  of  the  current 
reflect  essentially  number  density  and  electron  temperature 
fluctuations.'*^  The  natural  plasma  frequencies  were  obtained 
by  simultaneously  sweeping  the  triple  probe  and  the  Hall  probe 
through  the  plume  of  the  thruster.  The  triple  probe  was  aligned 
with  the  thruster  axis;  consequently,  the  uncertainty  of  the 
measurements  increased  because  it  was  no  longer  aligned  with 
the  ion  flow  vector  in  all  regions  of  the  plume.  The  uncertainty 
of  the  measurements  was  estimated  to  be-  ^^25%  for  T  and 
approximately  a  factor  of  2  for  n,  (Ref.  10). 

A  Nicoiet  320  digital  oscilloscope,  sampling  at  5  MHz,  re¬ 
corded  the  triple-probe  outputs,  resulting  in  the  maximum  ob¬ 
servable  frequency,  the  Nyquist  frequency^  of  2.5  MHz.  Of 
the  4000  data  points  obtained  by  the  digital  oscilloscope,  2048 
were  used  for  data  analysis  such  that  a  2"  data  point  fast 
Fourier  transform  (FFT)  algorithm  could  be  utilized.  Before 
computing  the  FFT,  the  time-averaged  value  of  the  signal  was 
subtracted  from  each  data  point  of  the  oscilloscope  trace,  and 
a  Hanning  window  function  was  applied  to  the  temporal  data 
to  reduce  leakage  effects.^^  Consistent  vdth  the  use  of  low-pass 
filters  to  minimize  aliasing  error,  the  frequency  response  of  the 
electronic  apparatus  was  tested  to  be  relatively  flat  up  to  ~1 
MHz.  All  power  spectra  shown  in  this  section  were  corrected 
for  the  frequency  response  of  the  instruments.  The  frequency 
response  of  the  triple  probe  is  expected  to  be  flat  up  to  the  ion 


plasma  frequency  (-100  MHz)  (Ref.  58),  which  was  well  be¬ 
yond  the  frequencies  investigated  (<2.5  MHz). 

Four  general  statements  can  be  made  about  the  measure¬ 
ments  taken  in  the  plume  of  the  10-kW  thruster.  First,  turbulent 
fluctuations  were  observed  only  at  high  values  of  Jj/m,  the 
minimum  value  ranged  from  40-60  kA^  s/g.  Second,  the  fluc¬ 
tuating  amplitudes  of  and  of  the  current  were  often  ob¬ 
served  to  be  a  significant  fraction  of  the  mean  values  (up  to 
50%).  Third,  the  power  spectra  of  were  essentially  identical 
to  that  of  the  current.  Finally,  the  power  spectrum  measure¬ 
ments  obtained  in  this  experiment  could  be  placed  into  two 
distinct  categories:  one  consistent  with  the  presence .  of  the 
GLHDI  and  the  other  consistent  with  the  presence  of  both  the 
GLHDI  and  the  ECDI. 

Before  showing  a  sample  of  the  first  type  of  power  spectra 
it  is  useful  to  discuss  where  in  frequency  space  the  peaks  are 
expected  to  occur.  To  assist  in  the  interpretation  of  the  spectra, 
the  linear  characteristics  of  the  microinstabilities  are  used 
throughout  this  section,  even  though  the  spectra  are  turbulent 
and  nonlinear  in  nature.  In  actuality,  only  the  peaks  in  the 
spectra  will  be  compared  with  linear  theory;  no  attempt  was 
made  to  describe  their  structure.  In  the  ion  rest  frame,  the 
frequency  at  which  the  power  spectrum  is  maximum  is  ex¬ 
pected  to  be  on  the  same  order  of  frequency  corresponding  to 
maximum  growth  (cf.  Table  1).  Because  the  ions  are  flowing 
with  respect  to  the  triple  probe,  the  peak  frequency  will  be 
Doppler  shifted.  The  following  equation  is  used  to  predict  the 
peak  frequency  in  the  laboratory  frame: 

U  =  \f*  +  [k^-W2v]\  (1) 

where  |  |  signifies  the  absolute  value,  Ui  is  the  ion  flow  vector, 
and /*  and  k*,  respectively,  represent  the  frequency  and  wave- 
vector  corresponding  to  maximum  growth  of  the  instability  in 
the  ion  rest  frame,  which  are  tabulated  in  Table  1. 

For  the  GLHDI,  the  maximum  temporal  growth  rate  in'  the 
ion  rest  frame,  occurs  at  a  frequency  near  the  lower  hybrid 
frequency,  but  the  magnitude  and  direction  of  k*  are  sensitive 
to  the  plasma  properties  and  are  not  known  accurately  for  the 
plasma  conditions  of  the  10-kW  thruster  plume.  Of  particular 
interest  is- an  estimate  of  the  largest  frequency for  which 
the  peak  in  the  turbulent  spectra  is  expected  to  be  observed  in 
the  laboratory  frame.  From  Eq.  (1),  it  is  seen  that  the  largest 
value  of  fub  occurs  for  the  largest  value  of  k*,  which  corre¬ 
sponds  to  k*ri^  ~  1  for  the  GLHDI  (cf.  Table  1).  Substimting 
iiito  Eq.  (1) 


/n.„-/LH  +  (H,/2777-^)  (2) 


where  rt,  =  Equation  (1)  can  be  written  as 


fmzx.  ^  _ 

/lh  V(2rcr/M,-) 


(3) 


Note  that  the  ratio  Ui/ilKTJMif^  is  typically  2-4  in  the  plume 
of  an  MPD  thruster  operating  on  argon. 

The  maximum  frequency  f,^  can  be  estimated  from  experi¬ 
mental  measurements  using  the  following  method:  /lh  and 
are  computed  from  the  measured  magnetic  field  and  T„  and 
the  ion  flow  velocity  is  estimated  from  the  following  relation': 


lh 


where  r„  and  are  the  anode  and  cathode  radius,  respectively, 
and  is  the  permeability  constant.  This  expression  for  the 
ion  flow  velocity  is  consistent  with  the  regime  of  thruster  op¬ 
eration  under  study  (i.e.,  when  the  electromagnetic  component 
of  thrust  is  much  greater  than  the  electrothermal  component), 
and  flow  velocity  measurements”  and  velocities  inferred  from 
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thrust  measurements’  in  the  lO-kW  thruster  also  agree  fairly 
well  with  the  previous  expression  (within  ~50%). 

A  sample  of  the  first  type  of  spectra  is  shown  in  Fig.  5  for 
argon.  The  normalized  power  spectrum  vs  frequency  is  plotted, 
where  the  normalized  power  spectrum  is  equal  to  Suif)/ 
Sii(/m),  and/„  is  the  frequency  at  which  the  power  spectrum 
Sii  is  maximum.  The  scatter  in  the  power  spectra,  due  to  the 
use  of  discrete  Fourier  transform  techniques,  is  not  physical. 
The  scatter  can  be  reduced  by  increasing  m,  which  is  chosen 
as  a  compromise  between  acceptable  statistical  scatter  and  fre¬ 
quency  resolution.'”'”  The  left-most  data  point,  at  '=5  kHz, 
represents  the  average  power  between  -100  and  <“10  kHz. 
Each  data  point  thereafter  represents  the  average  power  be¬ 
tween  ±5  kHz  of  the  corresponding  frequency.  In  addition, 
based  on  the  observation  of  oscillation-free  Hall  probe  mea¬ 
surements,  the  power  spectra  is  expected  to  fall  off  at  fre¬ 
quencies  less  than  on  the  order  of  100  Hz.  Indicated  on  the 
lower  frequency  scale  are  the  measured  natural  frequencies  of 
the  plasma,  and  along  the  top  of  the  plot,  associated  with 
the  GLHDI  and  for  the  ECDI  (to  be  discussed  later)  are 
'  identified.  The  essential  feature  of  these  results  is  the  presence 
of  large-amplitude,  naturally  occurring  turbulent  fluctuations 
near  the  lower  hybrid  frequency  (which  is  different  from  any 
of  the  other  natural  frequencies  by  at  least  two  orders  of  mag¬ 
nitude).  These  spectra  are  consistent  with  the  presence  of  the 
GLHDI  in  the  plume  of  the  10-kW  thruster.  Measurements 
•taken  in  the  plume  of  the  10-kW  thruster  operating  on  helium 
also  showed  this  same  phenomenon.'*’ 

Samples  of  Ihe  second  type  of  power  spectra  are  shown  in 
Figs.  6  and  7.  These  spectra  consist  of  the  same  broadband 
peak  near  the  Doppler-shifted  lower  hybrid  frequency,  but  with 
the  inclusion  of  high-power  spikes  at  integral  frequencies  from 
each  other.  These  spikes  generally  occur  at  much  higher  fre¬ 
quencies  than  the  Doppler-shifted  lower  hybrid  frequency  and 
are  indicative  of  the  ECDI’s  presence  in  the  plume  of  the  10- 
kW  thruster.  Although  difficult  to  view  on  a  logarithmic  scale, 
this  type  of  spectra  is  also  observed  in  Fig.  5,  with  faint  peaks 
at  1  and  2  MHz. 

Assuming  these  Spikes  are  a  manifestation  of  the  ECDI,  their 
location  in  frequency  space  can  be  described  by  the  following 
relation: 


Fig.  S  Comparison  of  the  turbulent  power  spectrum  to  measured 
natural  frequencies  of  the  plasma  and  to  the  frequencies  corre¬ 
sponding  to  the  maximum  linear  growth  rate  (argon,  J,  =  685  A, 
m  =  S  mg/s,  r  =  1.1  cm,  z  =  1.6  cm,  and  m  =  4). 


Fig.  6  Sample  turbulent  power  spectrum  of  the  second  type  ob¬ 
tained  in  the  plume  of  the  10-kW  MPD  thruster.  This  power  spec¬ 
trum  indicates  the  presence  of  the  ECDI  (argon,  predicted  fre¬ 
quencies;  -  1?0  kHz  and  -  700  kHz). 


[(UJV,)  -  1] 


(n  =  1,  2,  3,  . . .)  (5) 


where  "Fa  =  This  expression  ■was  obtained  by  sub¬ 

stituting/*  and  k*  from  Table  1  into  Eq.  (1),  and  using  the 
fact  that  k-Uc  >  0  is  required  for  the  instability  to  Operate. 
The  frequency  /,  can  be  predicted  from  the  measurements  of, 
B,  /,  ttj,  and  T^,  and  by  assuming  that  T,  =  T^,  estimating  the 
direction  of  U,  from  separate  flow  angle  measurements,"*’  and 
by  the  uSe  of  Eq.  (4)  to  calculate  t/,-.  When/,.]  was  calculated 
to  be  greater  than  the  Nyquist  frequency,  the  high-power  spikes 
were  not  seen,  and  when  it  was  predicted  that/.i  was  less 
than  the  Nyquist  frequency,  the  location  of  the  spikes  were 
observed  to  agree  with  Eq.  (5)  within  experimental  error 
(which  was  a  factor  of  2-3).  This  evidence  suggests  that  the 
ECDI  is  also  prominent  throughout  the  near-field  plume  of  the 
10-kW  thraster  (r  ^  3  cm,  z  <  D).  ’’ 

Finally,  it  should  be  stated  that  no  conclusions  can  be 
reached  about  microturbulence  occurring  near  the  other  natural 
frequencies,  because  they  were  not  in  the  viewing  range  of  the 
experiment.  Furthermore,  it  is  not  expected  that  this  turbu¬ 
lence,  if  it  exists,  will  be  consistently  Doppler  shifted  into  the 
viewing  range.  The  Doppler  shift  of  the  ion  cyclotron  har¬ 
monics  is  expected  to  be  too  small  to  be  seen  in  this  experi¬ 
ment,  because  typically  k*  is  much  smaller  than  the  inverse 
electron  gyroradius  (k*  «  r^)  when  the  GLHDI  is  also  ac¬ 
tive.®*  It  may  also  be  possible  for  microturbulence  occurring 
at  frequencies  near/p,  and  (in  the  ion  rest  frame) 

to  be  Doppler  shifted  into  the  viewing  range  of  the  experiment 


Frequency  (MHz) 


Fig.  7  Sample  turbulent  power  spectrum  of  the  second  type  ob¬ 
tained  in  the  plume  of  the  10-kW  MPD  thruster.  This  power  spec¬ 
trum  indicates  the  presence  of  the  ECDI  (argon,  predicted  fre¬ 
quencies:  ~  70  kHz  and/„-i  ~  830  kHz). 

if  k*-Uj<  0.  This  is  not  expected  to  be  observed  on  a  regular 
basis  because  for  the  microinstabilities  corresponding  to  these 
frequencies;  /,  «/*  and  /, «  |/c*  •  £/;i. 

V.  Measurement  of  the  Plasma  Dispersion  Relation 
Large-amplitude,  naturally  occurring  plasma  property  fluc¬ 
tuations  (i.e.,  plasma  turbulence)  are  the  nonlinear,  saturated 
state  manifestation  of  plasma  instabilities.  For  this  reason,  mi- 
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croinstabilities  in  the  plume  of  the  10-kW  thruster  were  iden¬ 
tified  by  comparing  only  the  gross  features  of  the  turbulent 
power  spectrum  with  linear  stability  theory.  A  more  direct  ap¬ 
proach  in  the  identification  of  microinstabilities  is  to  measure 
the  linear  dispersion  relation,  along  with  the  plasina  properties 
required  by  the  linear  stability  model,  and  compare  the  ex¬ 
perimental  results  with  theory. 

In  the  plume  of  the  10-kW  thruster,  where  »  1,  many 
practical  considerations  limit  the  scope  of  this  type  of  exper¬ 
iment.  Although  a  linear  stability  model  exists  for  the  MPD 
thruster  plasma  for  arbitrary  (Ref.  29),  that  theory  does  not 
include  gradient  effects  that  are  expected  to  be  important  at 
high  Pt-  Furthermore,  the  accuracy  of  plasma  dispersion  rela¬ 
tion  and  .plasma  property  measurements  are  limited  such  that 
a  careful  comparison  of  experiment  with  linear  theory  would 
not  be  fruitful;  indeed,  the  ion  temperature  and  ion  flow  ve¬ 
locity,  required  by  the  model,  were  not  measured  at  all. 
Throughout  this  study,  the  ion  temperature  was  assumed  to  be 
approximately  equal  to  I),  and  Ui  has  been  inferred  from  Eq, 
(4).  In  light  of  these  limitations,  it  was  desired  that  the  results 
of  this  experiment  would  provide  an  order-of-magnitude  type 
of  evidence  in  support  of  the  presence  of  the  GLKUDI  and 
ECDI  in  the  plume  of  the  10-kW  thruster. 

The  essence  of  this  experiment  consists  of  two  Langmuir 
probes  placed  a  finite  distance  apart.  These  probes  are  biased 
negative  with  respect  to  the  plasma  potential  to  measure  ion 
saturation  current  fluctuatioris,  which  are  assumed  to  reflect 
naturally  occurring  ion-number  density  fluctuations  only.  By 
measuring  the  phase  difference  between  the  fluctuations  at 
each  of  the  probes,  the  plasma  dispersion  relation  (i.e., 
along  the  line  of  the  probes  as  a  function  of  frequency)  can 
be  derived.  Technically,  a  nonlinear  theory  is  required  to  de¬ 
scribe  the  phase  information  obtained  from  this  experiment, 
but  as  mentioned  earlier,  it  is  commonly  observed  that  satu¬ 
rated-state  plasma  fluctuations  still  retain  much  of  their  linear 
characteristics.^^  Along  with  the  plasma  dispersion  relation 
measurements,  radial  profiles  of  the  plasma  properties  were 
measured  by  the  thmster-axis-aligned  triple  probe  and  the  Hall 
probe. 

Two  asymmetric  double  probes,  oriented  perpendicular  to 
the  thruster  axis,  were  used  to  measure  the  ion  saturation  cur¬ 
rent  fluctuations.''^  The  two  double  probes  were  placed  such 
.that  the  larger  electrodes  were  situated  nominally  5  mm  apart 
along  the  thruster  axis,  which  was  the  main  direction  of  the 
electron  drift  in  the  10-kW  thruster  plume.  Across  the  elec¬ 
trodes  of  the  two  probes  a  fixed  bias  (nominally  ~20  V)  en¬ 
sured  that  the  double  probes  operated  in  the  ion  saturation 
regime  and  ac  current  probes  were  used  to  record  the  ion  cur¬ 
rent  fluctuations.  Each  of  these  signals  was  then  amplified, 
filtered  (to  reduce  aliasing  error),  and  recorded  on  a  Nicolet 
320  digital  oscilloscope  sampling  at  a  2-MHz  rate.  The  data 
were  processed  in  the  same  manner  as  discussed  previously. 

Measurements  were  taken  at  various  regions  of  the  plume 
and  at  many  different  thruster  operating  conditions.  It  was  gen¬ 
erally  observed  that  the  phase  velocity  was  relatively  indepen¬ 
dent  of  frequency  over  the  frequency  range  of  100-500  kHz 
(0.7  ^  // /lh  ^  10),  which  was  the  region  where  the  mea¬ 
surements  were  most  accurate.  Typically,  ranged  from  =“8 
-15  km/s  for  50  ^  P,lm  (kA^  s/g)  ^  60,  and  from  '“20-30 
km/s  for  60  ^  /?//n  (kA^  s/g)  ^  90.  The  phase  velocity  was 
always  directed  in  the  +z  direction  for  regions  of  the  plume 
where  J.>0  and  <  0. 

A  typical  phase  velocity  measurement  is  shown  in  Fig.  8. 
The  error  associated  with  the  phase  velocity  measurement*’'  is 
due  to  statistical  scatter,  aliasing  error,  and  the  uncertainty  in 
the  phase  shift  due  to  the  electronics  (<5  deg  up  to  550  kHz). 
Although  the  signals  of  the  two  probes  were  very  coherent  at 
frequencies  less  than  500  kHz,  the  large  error  at  low  frequen- 
.  cies  was  due  to  the  large  fractional  uncertainty  of  the  fre¬ 
quency  and  phase  shift  measurement.  At  high  frequencies 
(>500  kHz),  statistical  scatter  and  aliasing  error  dominate.  The 


magnitude  of  and  the  observation  that  was  relatively 
constant  over  the  frequency  range  investigated  both  compare 
favorably  with  the  theoretical  characteristics  of  the  GLHDI. 
For  the  high-/8  region  of  the  MPD  thruster,  the  phase  velocities 
(in  the  ion  rest  frame)  are  expected  to  be  greater  than,  but  on 
the  order  of,  the  ion  thermal  velocity  (V,i  =  2.2[7)(eV')]‘'^  km/ 
s  for  argon)  and  relatively  independent  of  frequency.^"'''^  In  the 
laboratory  frame,  the  phase  velocity  will  be  increased  by  an 
increment  'Ui-kl\k\,  which  for  the  data  shown  in  Fig.  8  is 
approximately  12  km/s.  Assuming  7)  ~  T„  the  magnitude  of 
the  phase  velocity  is  definitely  near  the  theoretically  expected 
values. 

Many  times  the  power  spectra  of  the  ion  current  fluctuations, 
obtained  in  this  experiment,  exhibited  the  same  high-power 
spikes  seen  in  Figs.  6  and  7.  In  addition,  the  location  of  the 
spikes  in  the  frequency  domain  were  observed  to  agree  with 
Eq.  (5)  within  experimental  error.  By  investigating  the  mag¬ 
nitude  of  Vi  associated  with  these  spikes,  it  was  possible  to 
gain  further  evidence  for  the  existence  of  the  ECDI  in  the 
plume  of  the  10-kW  thruster.  In  theory,  the  phase  velocity  (in 
the  ion  rest  frame)  associated  with  each  harmonic  of  the  ECDI 
is  very  near  the  ion  thermal  velocity.  In  the  laboratory  frame, 
the  phase  velocity  is  expressed  as 

(6) 

This  expression,  which  represents  the  resonance  condition  be¬ 
tween  the  energy  absorbing  ions  and  the  negative-energy  Bern¬ 
stein  harmonics,^®  is  expected  to  be  very  resilient  among  the 
different  plasma  regimes  of  the  MPD  thruster.  In  fact,  in  its 


0  100  200  300  400  500 


Frequency  (kHz) 

Fig.  8  Typical  measurement  of  the  phase  velocity  in  the  plume 
of  the  10-kW  MPD  thruster  (argon). 
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Fig.  9  Typical  measurement  of  the  phase  velocity  in  the  plume 
of  the  10-kW  MPD  thruster  when  the  power  spectrum  exhibited 
sharp  peaks  similar  to  those  in  Figs.  6  and  7  (argon).  The  filled 
circles  correspond  to  the  peaks  in  the  power  spectrum. 
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turbulent  state,  as  in  this  experiment,  Eq;  (6)  is  also  expected 
to  be  valid.'*’ 

Figure  9  illustrates  a  typical  phase  velocity  measurement 
v/hen  the  high-power  spikes  were  present  in  the  power  and 
cross  spectra,  llie  filled  circles  at  the  frequencies  of  290  and 
580  kHz  represent  the  location  in  frequency  space  where  the 
spikes  were  observed.  The  phase  velocity  associated  with  the 
spikes  are  ‘=“9  km/s,  which  is  near  theoretical  expectations  (as¬ 
suming  Ti  ~  Tc).  Phase  velocity  measurements  in  agreement 
with  Eq.  (6)  7-12  km/s)  were  consistently  Observed  at 

the  frequencies  corresponding  to  the  spikes  in  the  power  and 
cross  spectra. 

In  summary,  order  of  magnitude  comparisons  of  linear  the¬ 
ory  with  phase  velocity  measurements  associated  with  the  nat¬ 
urally  occurring  plasma  fluctuations  were  'consistent  with  the 
presence  of  the  GLHDI  and  ECDI  in  the  plume  of  the  10-kW 
,MPD  thruster. 

VI.  Conclusions 

_  The  objective  of  this  research  was  to  provide  experimental 
evidence  for  the  existence  of  microinstabilities  in  a  10-kW- 
class  self-field  MPD  thruster.  The  lower-bound  of  the  range  of 
thruster  operating  currents  investigated  in  this  study  was  set 
by  the  requirement  that  the  electromagnetic  thrust  component 
be  much  greater  than  the  electrothermal  thrust  component;  the 
maximum  current  was  limited  by  the  so-called  onset  current. 

.  Because  of  the  lack  of  experimental  evidence,  the  theoretical 
characteristics  of  maihy  of  the  microinstabilities  that  may  exist 
in  the  self-field  MPD  thruster  were  reviewed.  It  was  concluded 
that  two  current-driven  microinstabilities  are  suspected  to  per¬ 
sist  throughout  the  plasma  of  the  self-field  MPD  thruster: 
GLHDI  and  ECDI.  From  the  comparison  of  the  location  of 
pealcs  in  the  power  spectra  (of  plasma  property  oscillations)  to 
the  linear  theory  of  the  GLHDI  and  ECDI,  experimental  evi¬ 
dence  for  their  existence  in  the  plume  of  the  lO-kW-class  self¬ 
field  MPD  thruster  was  obtained.  Additional  experimental  ev¬ 
idence  was  provided  by  the  measurement  of  phase  velocities 
associated  with  these  naturally  occurring  plasma  oscillations, 
which  were  also  observed  to  be  consistent,  within  experimental 
error,  with  the  linear  theory  of  the  GLHDI  and  ECDI. 
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